Spectroscopy, Interactions, and Level Splittings in Au Nano-particles. 
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We have measured the electronic energy spectra of nm-scale Au particles using a new tunneling 
spectroscopy configuration. The particle diameters ranged from 5nm to 9nm, and at low energies 
the spectrum is discrete, as expected by the electron-in-a-box model. The density of tunneling 
resonances increases rapidly with energy, and at higher energies the resonances overlap forming 
broad resonances. Near the Thouless energy, the broad resonances merge into a continuum. The 
tunneling resonances display Zeeman splitting in a magnetic field. Surprisingly, the g-factors 
(~ 0.3) of energy levels in Au nano-particles are much smaller than the g-factor (2.1) in bulk 
gold. 
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In nanoscale metallic particles, the electronic energy 
levels are quantized due to spatial confinement. Prior 
to the present work, this effect has been studied only by 
Ralph et al. 0], in small grains of Al. In this paper, we 
report tunneling data on energy spectra in Au nanopar- 
ticles, measured with a new tunneling configuration. We 
first show data which qualitatively confirm the results of 
Ralph et al. on the Al excitation spectra, but in Au parti- 
cles. Next, we confirm the role of the Thouless energy in 
the spectra of metallic particles. We also present new re- 
sults showing Zeeman splitting of energy levels analogous 
to the Zeeman effect in Al particles, but with g-factors 
much smaller than two. 

To be able to resolve the discrete energy levels in Au at 
~ 100 mK, the particle radius has to be less than approx- 
imately 10 nm. Fig. 1-a explains the steps in our sam- 
ple fabrication process. The first step is electron beam 
lithography. Using a PMMA bilayer resist technique, we 
define a resist bridge placed ~200nm above the Si wafer; 
this bridge acts as a mask. In step 2, we evaporate a 
15 nm thick film of Al, along the direction indicated by 
the arrow. Then we oxidize the surface of the film in 
50 mtorr of oxygen for 90 sec. In step 3, we deposit a 1 nm 
thick film of Au, along the same direction as in step 2. 
At this stage of film formation, Au forms isolated parti- 
cles, with a typical center-to-center spacing of 12 nm and 
a typical base diameter of 6 nm. Q Next we rotate the 
sample by ~90°. In step 4, we deposit the top electrode, 
along the direction indicated by the arrow. The top elec- 
trode is a bilayer, which has a ~3 nm thick layer of AI2O3 
at the bottom, and a 15 nm thick film of pure Al on top. 
The deposition angle is chosen so that the overlap be- 
tween the bottom and the top electrode is approximately 
20 nm. Typically, we grow many electrode pairs simulta- 
neously, and vary the overlap from to ^30 nm. Fig. 1-b 
is a schematic of the sample at a much larger magnifi- 
cation, showing that the particle is well screened by the 
electrodes. Because the oxide between the two electrodes 
is thick, current between the leads flows only through the 
particle(s), due to the exponential dependence of tunnel 



current on barrier thickness. 
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FIG. 1. a) Fabrication steps of our samples, b) Enlarged 
sketch of the particle, c) SEM image of a sample, top view. 

Fig. 1-c is an SEM image of a typical device. Au 
nanoparticles are visible on the top aluminum electrode 
and on the lower electrode far from the overlap region. In 
this work, we present data on particles with base diam- 
eter less than 10 nm, where their shape is roughly hemi- 
spherical. The area of the electrode overlap is such that 
in most cases, 0-3 particles are covered by the overlap. 
From the data, we can determine whether the tunneling 
current flows via only one particle, as explained in the 
next two paragraphs. The yield of single-particle sam- 
ples is about 30%. 

Fig. 2-a displays the I-V curve in sample 1 at 4.2 K. 
At this temperature, the energy levels are not resolved, 
and at first sight, the I-V curve is piece- wise linear. The 
electron transport at 4.2K is well described by the theory 
of single charge tunneling and the Coulomb staircase for 
tunneling via a single particle. j3J The junction capaci- 
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tances C\ and C%, and the background charge qo can be 
determined from the points where the I-V curve changes 
slope. From the theory ||, we can also estimate the 
junction resistances from the data. A summary of the 
parameters for this and two other samples are shown in 
Table 1. The capacitance per unit area in our junctions is 
« 50fF/fj,m 2 . Q We estimate the particle base diameter 
D and volume by assuming that the total capacitance of 
the hemispherical particle is equal to C\ + C2 ■ 
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FIG. 2. a) I-V curve of a typical device at 4.2K. b). 
I-V curve of the device at 30mK. c) First current step at 
three different refrigerator temperatures. The inset shows 
full-width-half-maximum vs. temperature of step a (after 
correcting for the capacitive division), d) Lineshape of the 
level when the leads are superconducting (H=0) compared to 
the lineshape when the leads are normal (H=1T). 

We apply a number of consistency checks to our mea- 
surements to ensure that we are measuring a single par- 
ticle. These consistency checks were introduced and ex- 
plained in the work by Ralph et al Q. Fig. 2-b presents 
the FV curve of the sample at 30 mK refrigerator tem- 
perature on an expanded scale. At low voltages, current 
increases in discrete steps. Fig. 2-c shows step a in detail 
at three different refrigerator temperatures. The curves 
between the points are fits to the Fermi distribution. The 
strong temperature dependence of the line-shape demon- 
strates that tunneling is via a discrete quantum level in 
the particle. After correcting for the capacitive division 
of voltage, the full-width-half-maximum (FWHM) of the 
peak is linear with the refrigerator temperature, as shown 
in Fig. 2-c, with the slope of 3.8fcsT, which is close to 
the expected slope of 3.5fcsT. At very low refrigerator 
temperatures, FWHM is temperature independent, and 
the base electron temperature is ^70 mK. A magnetic 
field of IT is applied to supress superconductivity in the 
leads. 

In Fig. 2-d we show the line-shape of the level when 
the Al leads are superconducting and when the supercon- 
ductivity is suppressed with a field of IT applied parallel 



to the film plane. In zero field, the line-shape traces the 
BCS density of states in the aluminum leads. The shift 
of the level which occurs when the leads change from su- 
perconducting to normal is proportional to the BCS gap 
in Al, and it depends on the ratio of the capacitances 
between the particle and the two leads. The shifts of the 
levels which we observe when Al is driven normal (not 
shown) indicate that all the levels have the same capaci- 
tances to the leads, as expected if all resonances occur in 
a single grain. 

Fig. 3 shows the excitation spectra (dl/dV) in three 
samples, in a field of IT to suppress the superconduct- 
ing gap in the leads. Peaks occur at the energy levels 
of the particle. For each of the samples, the corrective 
factor between voltage and particle energy is shown in 
the figure. In sample 2, not all the levels correspond to 
the same number of electrons on the particle. 

In the single-electron model, energy levels near the 
Fermi energy should be distributed with only small fluc- 
tuations between the successive level spacings (spectral 
rigidity). The spectra shown in Fig. 3 are very differ- 
ent from such single-electron excitation spectra. At low 
energies, spectroscopic peaks are sharp and well sepa- 
rated. As the excitation energy is increased, the density 
of tunneling resonances increases, and the peaks have 
a tendency to cluster. Eventually, we can not resolve 
individual energy levels anymore; however we can still 
resolve envelopes (broad resonances) of many overlaping 
tunneling resonances. The broad resonances increase in 
width as the energy increases, and finally, broad reso- 
nances merge forming a continuum in the tunneling den- 
sity of electronic states. 

The increase in the density of tunneling resonances 
with energy seen in Fig. 3 can be explained by electron- 
electron interactions. In the non-equilibrium model due 
to Agam et al. ^ , the same single electron state can oc- 
cur with different excited configurations of the other elec- 
trons in the particle. The theory of Altshuler et al. Q] 
describes the derealization (in Fock space) of a quasipar- 
ticle by the creation of electron-hole pairs. This leads to 
an increase in the density of discrete levels with increas- 
ing energy and to the existence of quasiparticles with fi- 
nite lifetime above a certain energy. The highest energy 
at which the quasiparticles can be resolved is predicted 
to be the Thouless energy. ]8| 

The idea that the Thouless energy sets the limit of the 
observability of the discrete energy spectrum in an in- 
teracting electron system has been first experimentally 
studied by Sivan et al. Q , on a disordered semiconduct- 
ing quantum dot. In that work, the authors' data anal- 
ysis shows that there has to be a finite energy above 
which the energy levels in a quantum dot can not be re- 
solved. In our particles, unlike semiconducting quantum 
dots, the charging energy is so much larger than the level 
spacing that we can, for the first time, access and display 
the entire progression, out to the high energy horizon, all 
within the first step of Coulomb blockade. 
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FIG. 3. Excitation spectra in three different samples, at 
T — 30mK and H = IT. In order to convert from voltages 
to particle eigenenergies, voltage must be multiplied by the 
appropriate correction factor, indicated for each graph. 

Because the density of observed peaks increases rapidly 
with energy, we rely on spectral rigidity and estimate the 
single-electron level spacing 8 as the distance between the 
two lowest energy peaks. The agreement between the 
measured level spacing and the calculated one (from the 
grain volume inferred from the capacitance) is within a 
factor of 2, as shown in Table 1. We note that in sample 3, 
the distance between the low energy clusters is approxi- 
mately the single-electron level spacing, consistent with 
the theories of electron-electron interactions. 

If we assume that the electronic motion in the particle 
is ballistic, and that the surface scattering is diffusive, 
then the mean-free path I of an electron confined in the 
particle is approximately D. To diffuse through the en- 
tire particle volume however, an electron should scatter 
from the particle surface about three times. Thus, we 
estimate the electronic traversal time as 3D/vp (essen- 
tially, the factor of 3 is the number of dimensions). The 
Thouless energy is given approximately as Et = Tivp /3-D 
= 37meV, 75meV and 40meV, for samples 1,2, and 3, re- 
spectively. It can be seen from samples 1 and 2 that these 
estimates of the Thouless energy are consistent within a 
factor of 2 of the bias voltage at which the spectra be- 
come independent of energy. 

Note that the agreement between the theory and our 
data is better if we compare the bias voltage with the 
Thouless energy, i. e., if we do not correct the Thouless 
energy with capacitive division prefactors. Our expla- 
nation is that the highly excited states predicted by the 
Agam model can initiate the formation of electron-hole 
pairs. The maximum number of a particle's excited states 
which can form when electrons exit the particle, is given 
by the bias voltage divided by the level spacing. This 
number is ~10, 12, and 20 for samples 1, 2, and 3, re- 
spectively. H Thus, the role of large bias voltage (due to 
large charging energy) is to reduce the range of energies 
where the levels can be resolved. If qo ~ e/2, so that 
there is little Coulomb blockade, we would expect to see 



many more resolved levels, as found by Ralph et al. in 
gated aluminum particles. 

We have measured the evolution of the spectrum with 
applied magnetic field. In a metallic particle, every or- 
bital state has a 2-fold spin degeneracy. Let n be the 
number of electrons on the particle after an electron has 
tunneled onto (or off) the particle. If n is even, the total 
spin in the ground state is zero in the simplest model of 
weakly interacting electrons, and the first tunneling res- 
onance shifts but does not split with magnetic field. If n 
is odd, there is an unpaired spin, and all the levels split 
with magnetic field. (Because the particle is small, the 
effect of the field on electronic orbits is negligible.) 

Zeeman splittings in agreement with this simple model 
have been measured in Al nanoparticless [jj] and some 
carbon nanotube ropes JlO| |. However, in ferromagnetic 
particles, and, surprisingly, also in large semiconduct- 
ing quantum dots and some carbon nanotubes, the ex- 
pected spin degeneracy is not observed experimentally. 
p"T[ In this work, we do measure Zeeman splitting in 
Au nanoparticles, confirming that spin degeneracy is 
commonly observed in simple metals, as expected, even 
though it is not observed in some other systems. 

Fig. 4 is the gray-scale image of the magnetic field de- 
pendence of the energy spectrum in sample 1, at negative 
bias voltage, where the particle is further away from equi- 
librium than for positive bias. The low energy levels split 
with magnetic field, and the splitting is linear with field. 
Because the lowest level splits into two, we conclude that 
the number of electrons is odd. From the data, we obtain 
that the g- factor in this resonance is g — 0.28. Sample 3 
has similar splittings, with a g-value of 0.44. In sample 2, 
the lowest two peaks do not split implying that the num- 
ber of electrons on the particle is even. The splitting of 
the higher energy peaks in sample 2 could not be resolved. 

A <?-factor of ~ 0.3 in Au nanoparticle eigenenergies is 
much smaller than the g- factor of 2.1 determined by elec- 
tron spin resonance measurements in bulk Au. |l2j ] We 
explain the weak Zeeman splitting of Au nanoparticle 
eigenenergies by the strong spin-orbit (s-o) interaction in 
Au. The explanation is corroborated by the work by Sali- 
nas et al. |ll| , where eigenenergies of superconducting Al 
nanoparticles doped with 4% of Au had g- factors in range 
0.5-0.8, instead of 2, which is the g-factor in pure Al. The 
effect of s-o interaction on the single-electron eigenstates 
JLlj is to mix up-spin and down-spin polarizations, which 
leads to reduced g- factors. 

The discrepancy between the g-factors determined by 
spin resonance experiments in bulk Au ]l2| and the g- 
factors of ~ 0.3 obtained here by tunneling spectroscopy 
implies that for a high-Z element such as Au, these two 
techniques measure two different quantities, which can 
not be directly compared. Further theoretical study 
is needed to explain this discrepancy and derive the 
g-factors in Au nanoparticles. We speculate that the 
nanoparticle geometry (very large surface to volume ra- 
tio) plays an essential role in g-factors reduction, pres- 
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sumably because the spin-flip scattering in Au nano- factors in bulk Au 
particles occurs at the particle surface. 

Note that in Fig. 4, resonance (3 does not split with 
magnetic held, and resonances 7 and 7' split with mag- 
netic held and cross; the analysis of these effects is beyond 
the scope of the present paper. |1| 
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FIG. 4. Magnetic field dependence of the energy spectrum 
in sample 1. Top schematic is a guide to the eye. 

In conclusion, we have performed spectroscopy of dis- 
crete energy levels in Au nanoparticles, with a new tun- 
neling configuration. Data for some of our samples in- 
clude the entire spectrum starting from single electron 
levels at low energy, progressing to a regime where only 
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be resolved, and ending in a continuum. We have ob- 
served Zeeman splitting under applied magnetic held and 
found very small g-factors, much smaller than the g- 
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TABLE I. Ci, C2, Ri, and R2: junction capacitances and resistances determined from Coulomb staircase at 4.2K. D: particle 
base diameter estimated from Ci + C2 assuming hemispherical shape. 5: Estimated spacing beween electron-in-a-box levels, 
based on particle volume. 5: measured level spacing. Et'- The Thouless energy of particles, estimated as ?ivf/SD', where D' 
is the particle base diameter corresponding to measured level spacing 8. g: determined from Zeeman splitting. 
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